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ABSTRACT 



We present new observations made with the IRAM 30m telescope of the J=l-0 and 3-2 lines of HCN and HCO^ used to probe 
the dense molecular gas content in a sample of 17 local luminous and ultraluminous infrared galaxies (LIRGs and ULIRGs). These 
observations have allowed us to derive an updated version of the power law describing the correlation between the FIR luminosity 
(Lfir) and the HCN(I-O) luminosity (^hcnci-o)) °f local and high-redshift galaxies. We present the first clear observational evidence 
that the star formation efficiency of the dense gas (SFEjense). measured as the ^fir/^hcn(i-o) r^tiO' is significantly higher in LIRGs and 
ULIRGs than in normal galaxies, a result that has also been found recently in high-redshift galaxies. This may imply a statistically 
significant turn upward in the Kennicutt-Schmidt law derived for the dense gas at Lfir ^ 10" Lq. We have used a one-phase Large 
Velocity Gradient (LVG) radiative transfer code to fit the three independent line ratios derived from our observations. The results 
of this analysis indicate that the [HCN]/[HCO^] abundance ratios could be up to one order of magnitude higher than normal in a 
significant number of LIRGs and ULIRGs of our sample. An overabundance of HCN at high Lfir implies that the reported trend in the 
ipiR/^HCN ratio as a function of Lfir would be underestimating a potentially more dramatic change of the SFEdcnsc- Results obtained 
with two-phase LVG models corroborate that the LJ_ 



HCN(l-O) 



■to-Mjensc conversion factor must be lowered at high Lfir. We discuss the 
implications of these findings for the use of HCN as a tracer of the dense molecular gas in local and high-redshift luminous infrared 
galaxies. 

Key words, galaxies: evolution - galaxies: ISM - galaxies: starburst ~ infrared: galaxies - ISM: molecules - radio lines: galaxies 
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The question of how the star formation rate (SFR) in galaxies 
scales with the density of the star-form ing gas is a key prob- 
lem in extragalactic research. I Schmidt! ([1959) first postulated 
that the SFR per unit volume (psfr) should vary as a power 
law of index of the gas volume density (pgas): Psfr Pg-^s- 
The translation of the Schmidt law in terms of the coiTespond- 
ing surface densities of SFR (Esfr) and gas (Egas) would imply 
SsFR oc S^jij, if we assume a roughly constant scale-height of 
the star-forming gas in galaxies. To observationally validate the 
Schmidt law requires in practice finding fair quantitative tracers 
of EsFR and Eg^s. Above a certain density threshold, determined 
by large-scale gr avitational instability of disks (e.g., iKennicutj 
ll989l) . lKennicuttl ( il998,) found that the disk-averaged SFRs and 
gas densities of a sample of ~100 galaxies were well repre- 
sented by a Schmidt law wit h an index ~ 1.4, the so called 
Kennicutt-Schmidt (KS) law . iKeniiicuttl (flgH) used CO and HI 
data to derive Sgn^. lGao & S olomon (2004a, b, hereafter GS04a 
and GS04b), using a sample of 65 galaxies, derived a similar su- 
perlinear correlation between the infrared and the CO(l-O) line 
luminosities, taken as proxies for the SFR and the total molec- 



ular gas content, respectively: Ljr oc L' 



1.4-1.7 
CO(l-O)' 



Interestingly, 



they found instead a tight linear correlation (i.e., ~ 1) over 



* Based on observations carried out with the IRAM 30-meter tele- 
scope. IRAM is supported by INSU/CNRS (France), MPG (Germany) 
and IGN (Spain). 



3 decades in Lir between Lir and the luminosity of the HCN(1- 
0) line i'HCN(i-o)' latter being a tracer of dense molecular 
gas (hht > 10"^ cm"^). More recentlv ^Vu et al.l (|2005|) have ex- 
tended a similar correlation to much smaller scales by observing 
Galactic dense cores in HCN(l-O) emission. 

On theoretical grounds it is expected that the power law of 
the KS law should be close to 1.5. If star formation proceeds due 
to small-scale gravitational collapse of the gas, an index = 1 .5 
comes naturally, provided that a constant fraction of the gas 
forms stars in a free-fall time (e . g . . iLarsonll 1 9 8 8 ) . Alternatively, 
in a Toomre-stable galaxy disk, an index A^ = 1 .5 is expected if 
a constant fraction of the disk gas forms stars per unit dynamical 
time-scale, this being determined by the rotation period (e.g., 
Elmegreen 2002). More recendy, Krumholz & McKee (2005|) 
have advanced a model for star formation regulated by super- 
sonic turbulence where psFR would be a power law of the mean 
density of the gas, p^^^, with an index A^ ~ 1.5. 

Ba sed on the model developed by Krumholz & McKei 
(I2005L see also disc u ssion in .Krumholz & Tan .200' 
iKrumholz & Thompson! (!2007!) have provided a frame- 
work describing to what extent the derived correlation between 
Psfr and Pg^j, would depend on the molecular line tracer used 
as a proxy for p„^^. Line transitions of high effective critical 
densities (weff ~ 10"*"^ cm"-*; e.g, like the low-J rotational lines 
of HCN and HCO^) would only trace high-density peaks, and 
they would be thus insensitive to the bulk of the molecular 
gas in galaxies. On the contrary, transitions characterized by 
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Table 1. Main properties and IRAM 30m telescope observational results of our sample of LIRGs and U LIRGs. 0co is the molecula r 
source size at FWHM derived from interferometric CO observations found in the literature (mainly from' Dovynes & S olomonl l998l) . 
Velocity int egrated line intensities are g iven in antenna temperature scale (T*). 1-cr uncertainties are derived from the basehne fits. 
Data from lGracia-Carpio et all (l20od) . 



Source 


D A 

K.A. 


Decl. 




Zco 


n 


'HCN(l-O) 


'HCN(3-2) 


T 

JHCO+(1-0) 


T 

JHC0+(3-2) 




(J zOUO) 


(J 2000) 


(Mpc) 




( ) 


(Kkms ) 


(Kkms ) 


..—1 \ 

(Kkms ) 


(Kkms ) 


IRAS 17208-0014 


17 23 21.90 


-00 17 00.1 


187 


0.04288 


1.70 


2.19 ± 0.16" 


4.73 ± 0.43 


1.49 ± 0.15" 


3.06 ± 0.33" 


Mrk231 


12 56 14.20 


56 52 25.9 


186 


0.04220 


0.85 


1.39 ±0.08 


1.95 ±0.22 


1.34 ±0.12" 


2.04 ± 0.25" 


IRAS 12112+0305 


12 13 46.00 


02 48 41.0 


331 


0.07291 


1.00 


0.62 ±0.12" 


< 1.70 


0.37 ± 0.09" 




Arp 220 


15 34 57.20 


23 30 11.5 


80 


0.01818 


1.80 


8.16 ±0.17 


18.04 ±0.51 


3.77 ±0.21" 


4.57 ± 0.23 


Mrk 273 


13 44 42.10 


55 53 13.1 


166 


0.03776 


0.73 


1.11 ±0.13 


3.07 ±0.61 


1.15 ±0.12" 


2.90 ± 0.23" 


IRAS 23365+3604 


23 39 01.30 


36 21 10.4 


280 


0.06438 


0.95 


0.39 ± 0.07" 


<0.87 


0.26 ± 0.06" 




UGC 05101 


09 35 51.60 


61 21 11.6 


173 


0.03931 


3.50 


1.40 ± 0.14 




0.75 ±0.13 


2.13 ± 0.24 


VII Zw 31 


05 16 46.70 


79 40 12.0 


238 


0.05429 


2.24 


0.48 ± 0.07 




0.56 ± 0.08" 


1.06 ±0.26 


NGC 6240 


16 52 58.80 


02 24 03.8 


106 


0.02448 


2.00 


2.18 ± 0.17 


9.18 ±0.68 


4.06 ± 0.21" 


4.24 ± 0.52" 


Arp 55 


09 15 55.20 


44 19 54.7 


176 


0.03984 


4.50 


0.74 ± 0.07" 




0.69 ± 0.09" 


0.95 ± 0.23 


Arp 193 


13 20 35.30 


34 08 24.6 


104 


0.02335 


1.50 


1.04 ± 0.09 


1.13 ±0.21 


1.67 ±0.13" 


3.11 ±0.22" 


NGC 695 


01 51 14.30 


22 34 56.2 


136 


0.03245 


4.00 


0.43 ± 0.08" 




0.62 ± 0.08" 


<0.69 


Arp 299 A 


11 28 33.50 


58 33 45.3 


47 


0.01044 


5.00 


2.04 ± 0.11" 


1.58 ±0.38 


3.88 ±0.16" 


4.72 ± 0.44" 


Arp299B+C 


11 28 30.80 


58 33 48.3 


47 


0.01044 


7.00 


1.29 ± 0.09" 


<0.72 


2.09 ±0.16" 


1.64 ±0.24" 


NGC 7469 


23 03 15.60 


08 52 26.3 


65 


0.01643 


4.24 


1.85 ± 0.09" 


2.76 ± 0.34 


2.70 ±0.10" 


2.08 ± 0.42 


Mrk 331 


23 51 26.80 


20 35 10.0 


72 


0.01805 


4.00 


1.35 ±0.10" 


1.16 ±0.27 


1.74 ± 0.09" 


2.18 ±0.41" 


NGC 7771 


23 51 24.90 


20 06 42.6 


56 


0.01428 


4.00 


3.81 ± 0.13" 


3.33 ±0.58 


3.67 ± 0.30" 


1.98 ±0.37 



lower effective critical densities (Weff ~ lO^^-'cm"^; e.g., like 
the low-J rotational lines of CO) would be better tracers of 
the median density of the bulk of the gas. In this scenario, the 
different power indexes of SFR so far obtained using either CO 
or HCN could be explained without violating the universality 
of the SFR law in ga l axies described above. In particular, 
iKrumholz & Thompson! (l2007h argu e that the linear correlation 
between Ljr and iHCN(i-o) found by GS04b cannot be taken as 
firm evidence that the objects detected by the HCN(l-O) line 
represent a physically distinct star forming unit in molecular 
gas. 

The use of HCN lines as a true quantitative tracer of the 
dense molecular gas in galaxies can face difficulties if the ex- 
citation conditions and/or the chemical environment of molec- 
ular gas depart from normalcy. This can significantly change 
the conversion factor between the luminosity of HCN lines and 
the mass of dense molecular gas. From the observational point 
of view, there is mounting evidence that HCN lines can be 
overluminous with respect to the lines of other dense molec- 
ular gas tracers, like HCO^, in the cir cumnuclear disks of 
Seyferts jKohno et all l200lt lUsero et all 12004; Kohno 2005). 
A significant percentage of luminous and ultraluminous in- 
frared galaxies (LIRGs and ULIRGs) has also been reported 
to show high HCN/HCO^ intens ity ratios dGracia-Carpio et al] 
|2006t llmanishi et al."2006', '2007^. A similar result has been re- 
cently found by Garcia-Burillo et al. (2006) in the BAL quasar 
APM 08279+5255 at z ~ 4. The origin of overluminous HCN 
lines, in the terms described above, is still unclear and sev- 
eral theoretical scenarios have been advanced in the litera- 
ture. X-rays may significantly enhance the abundance of HCN 
relative to other molecular species like HCO^ in enshrouded 
AGNs, where X-ray D ominated Regions (XDR) can develop 
( Lepp & Dalg arnd 4996). More r ecently, iMeiierink & SpaansI 
(l2005h andlMeiierink et al. (|2007|) have proposed that Photon 
Dominated Regions (PDRs) are more efficient than XDRs in 
elevating HCN-to-HCO^ ratios in starburst galaxies. Hot core- 
like chemistry in starbursts have also bee n invoked as a mech - 
anism responsible of enhancing HCN dLintott & Vitil l2006h . 



Furthermore, instead of being collisionally excited, HCN lines 
might reflect the pumping by IR photons. These conditions can 
prevail in the molecular circumnuclear disks around the strong 
mid infrared s ources typically found in AGNs (Aalto et al. 1995 ; 
Garcia-Burill o et alJ l2006l: iGuelin et al.ll2007t IWeiB et alJl2007 : 
Aalto et aL,2007bllal) . 

The caveats on the use of HCN as the only standard tracer 
of the dense molecular gas in galaxies call for the observation of 
other molecular species and transitions, in particular in LIRGs 
and ULIRGs. This question is central to disentangling the dif- 
ferent power sources of the huge infrared luminosities of these 
galaxies. In this paper we present new observations made with 
the IRAM 30m telescope of the J=l-0 and 3-2 lines of HCN 
and HCO^ used to probe the dense molecular gas content of a 
sample of 17 LIRGs an d ULIRGs. Preliminary re sults of this 
work were published by iGracia-Carpio et al.l (|2007|) . These ob- 
servations, t hat complement t he first HCO^ survey published 
by iGracia-Carpio et all (l2006l hereafter GC06) of LIRGs and 
ULIRGs, are used to derive a new version of the power law de- 
scribing the correlation between LpiR and i[jcN(i_o) from normal 
galaxies (Lfir < 10" Lq) to high-z galaxies. We present the first 
clear observational evidence that the star formation efficiency of 
the dense gas, measured as the L^iiYil L'^^^My-i)) ratio, is signifi- 
cantly higher in LIRGs and ULIRGs than in normal galaxies. 
We also find that [HCN]/[HCO^] abundance ratios could be up 
to one order of magnitude higher than normal in a significant 
number of LIRGs and ULIRGs. We discuss the implications of 
these findings for the use of HCN as a tracer of the dense molec- 
ular gas in local and high-redshift luminous infrared galaxies. 

2. Observations 

The new HCN(l-O), HCO+(1-0), HCN(3-2) and HCO+(3-2) 
observations were carried out in five observing runs between 
December 2005 and November 2006 with the IRAM 30m tele- 
scope at Pico de Veleta (Spain). The sample consists of 17 
LIRGs and ULIRGs selected to cover homogeneously the Ljr 
range between 10"^ Lq and 10'^'^ Lq. All galaxies are located 
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Fig. 1. HCN, HCO+ and CO spectra observed with the IRAM 30m telescope. The HCN(3-2), HCO+(3- 
sities have been scaled by the factor indicated in the panel. The velocity windows used to calculate the 
are highlighted in grey. Galaxies are ordered with decreasing Lfir- 



2) and CO(l-O) line inten- 
line areas given in Table [1] 



at distances larger than 50Mpc to be confident that the total 
emission of the molecular gas can be measured in a single point- 
ing: FWHM(90GHz) ~ 28" = 7kpc and FWHM(260 GHz) 
~ 9" = 2.5 kpc, at 50 Mpc. The 3 mm and 1 mm SIS receivers of 
the 30m telescope were tuned to the redshifted frequencies of the 
lines. The velocity range covered was 1800 kms"' for the 3 mm 
lines and 1200 km s"' for the 1 mm lines. The wobbler switching 
mode was used to obtain flat baselines, that was almost always 
the case for the 3 mm observations. Only those individual spectra 
at 1 mm that showed a flat profile, where an order zero polyno- 
mial (i.e., a constant) had to be subtracted, have been retained; 
others showing a tilted profile were rejected and have not been 



included in the resulting average spectrum. The velocity ranges 
used to fit the flat baselines were chosen to be identical to those 
used for the higher signal-to-noise CO profiles obtained as part 
of this survey for all the galaxies of our sample. 

Typical system temperatures during the observations were 
~130K at 3 mm and ~700K at 1 mm. All receivers were used in 
single side-band mode (SSB), with a high rejection of the image 
band: >12dB at 1mm and >20dB at 3mm. The latter assures 
that the calibration accuracy for the bulk of our data is better 
than 20%. Pointing of the 30m telescope was regularly checked 
every 2 hours by observing nearby continuum sources; we found 
an average rms pointing error of 2"-3" during the observations. 
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Fig. 1. Continued 



When occasionally a larger poinling error was found during an 
observation, some of the spectra previously taken were dropped, 
the pointing was corrected, and additional spectra were taken 
towards the source. 

Throughout the paper, velocity-integrated line intensities (/) 
are given in antenna temperature scale, T*. The T* scale re- 



lates to the main beam temperature scale, Tmb, by the equa- 
tion r^nb - (Feff/Bgff)T*, where Feff and Beif ^ire the forward 
and beam efficiencies of the telescope at a given frequency. For 
the IRAM 30m telescope i^eff/Beff = 1-22 (1.96) at 86 GHz 
(260 GHz) and 5/r„,b = 4.95 Jy K'' . The velocity window used 
to calculate /hco+ and /hcn in the two rotational lines has been 
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Fig. 1. Continued 



defined on a case-by-case basis by using the CO(l-O) line pro- by Ho - 71 kms ' Mpc ' and Q.^ - 0.27 dS pergel et al.ll2003l) . 

files. Molecular line luminosities (L') were computed in units Results are summarized in Tables [1] and |2l All the HCN, HCO^ 

of L' = Kkm s ' pc^ according to Equation 1 of lGS04a. where and CO line profiles used in this work are displayed in Fig.[T] 
the K scale here corresponds to the brightness temperature av- 
eraged over the size of the source. Luminosity distances have 
been derived assuming a flat A-dominated cosmology described 
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Table 2. Derived luminosities and surface densities of our sample of LIRGs and ULIRGs. L' - Kkm s ' pc^. ^ Data from lGC06l 
corrected to the different cosmology adopted in this paper and to the finite source size of the molecular gas estimated from interfer- 
ometric CO observations found in the literature (mainly from Downes & Solomon 1998). 



Source 


^FIR 

(10" Lo) 


LiR 

(10" Lo) 


-^HCN(l-O) 

(10** L') 


■^HCN(3-2) 

(10** L') 


^HCO+(l-0) 

(lO** L') 


^HC0+(3-2) 

(10** L') 


SsFR 

(Moyr"' kpc"^) 


^dense 

(Mopc-2) 


IRAS 17208-0014 


18.5 


25.6 


18.4 


± 1.3" 


7.2 ± 0.7 


12.4 


± 1.2" 


4.6 ± 0.5" 


235 


10400 


Mrk 23 1 


17.3 


35.0 


11.6 


±0.7 


2.9 ± 0.3 


11.0 


± 1.0" 


3.0 ± 0.4" 


880 


26400 


Tn AO 1 1 1 ■ none 

IRAS 12112+0305 


15.1 


21.9 


15.8 


± 3.1" 


< 7.5 


9.3 


± 2.3" 




220 


10400 


Arp 220 


11.7 


15.5 


12.7 


±0.3 


5.1 ±0.1 


5.8 


±0.3" 


1.28 ± 0.06 


600 


^ n'^ nn 

29200 


Mrk 273 


9.9 


14.9 


7.4 


±0.9 


3.6 ± 0.7 


7.6 


±0.8" 


3.4 ± 0.3" 


620 


'^nonn 

20800 


IRAS 23365+3604 


9.7 


15.2 


7.2 


± 1.3" 


< 2.8 


4.8 


± 1.1" 




190 


6340 


UGC 05101 


7.7 


9.9 


10.2 


± 1.0 




5.4 


±0.9 


3.0 ± 0.3 


18 


1070 


VII Zw 31 


7.4 


9.7 


6.5 


± 1.0 




7.5 


± 1.1" 


2.6 ±0.6 


34 


1330 


NGC 6240 


4.6 


7.1 


6.0 


±0.5 


4.7 ±0.4 


11.0 


±0.6" 


2.2 ±0.3" 


140 


7960 


Arp 55 


4.1 


5.4 


5.7 


±0.5" 




5.2 


± 0.7" 


1.5 ±0.4 


7 


460 


Arp 193 


3.5 


4.8 


2.7 


±0.2 


0.6 ±0.1 


4.4 


±0.3" 


1.5 ±0.1" 


56 


1930 


NGC 695 


3.4 


4.6 


2.0 


±0.4" 




2.8 


±0.4" 


<0.7 


16 


410 


Arp 299 A 


2.6 


4.0 


1.14 + 0.06" 


0.20 ± 0.05 


2.15 


± 0.09" 


0.60 ± 0.06" 


60 


1190 


Arp299B+C 


1.7 


3.1 


0.75 ± 0.05" 


< 0.11 


1.19 


± 0.09" 


0.25 ± 0.04" 


20 


400 


NGC 7469 


2.3 


3.7 


2.0 


±0.1" 


0.62 ± 0.08 


2.9 


±0.1" 


0.46 ± 0.09 


33 


1270 


Mrk 331 


1.7 


2.6 


1.8 


±0.1" 


0.32 ± 0.07 


2.2 


±0.1" 


0.6 ±0.1" 


27 


1220 


NGC 7771 


1.6 


2.0 


3.0 


±0.1" 


0.54 ± 0.09 


2.9 


±0.2" 


0.33 ± 0.06 


7 


545 



3. A new HCN/HCO+ survey in LIRGs/ULIRGs: 
comparison with previous data 

The new data presented in this work represent a significant im- 
provement in the completeness and quality of the available sur- 
veys of the dense molecular gas cont ent of L IRGs and ULIRGs. 
Contrary to the previous survey of GS04a', based on a single 
transition of HCN (the J=l-0 line), the present data set includes 
another tracer of the dense molecular gas, HCO^, and expands 
the number of observed transitions to two per species. More 
relevant to the discussion of this paper, and as argued below, 
the new HCN(l-O) data partly change the results presented in 
lGS04bl concerning the constancy of the star formation efliciency 
in galaxies. 

Fig. |2h represents the Lir/L^^j^i^j luminosity ratio as a 
function of Lir derived from the data published by ,GS04b. for 
a sample of 65 targets including normal galaxies, LIRGs and 
ULIRGs (hereafter Lir refers to LiR(8-1000yum)). This plot is 
virtually identical to Fig. 2a of GS04b; only a few data points 
have moved in the diagram due to the use of more recently de- 
te rmined v alues for the distances to the galaxies. The main result 
of lGS04bf s paper is echoed by Fig. |2^: there is no strong trend 
of the Lir/LJ^(-,j^^j gj luminosity ratio with Ljr. A power law fit 



to these data gives Lir/L^(,^^j „ 



oc L, 



0.11-0.16 



, depending on 



whether limits on i^cN(i-O) included (a = 0.11 +0.05) or ex- 
cluded {a - 0. 16 +0.05) in the fit. In terms of average values, the 
i'lR/^HCNCi-O) luminosity ratio changes from ~900LqL'"' for 
normal galaxies (Lir < lO" L©) to -IIOOLqL' ' for ULIRGs 
(Lir > 8 X 10" L©), i.e., a marginally significant ~20% in- 
creas43. 



Note that the value of Lir/LJ^j^^j^j = IAQLqL' ' derived by 



lGS04 b for normal galaxies, and quoted in Table 2 of their paper, is in- 
consistent with the data values listed in Table 1 of GS04b; based on 
these data we rather derive Lir/LJj^j^,,, ~ 900 LoL'"' for normal 
galaxies. To ease the comparison with iGS04b results, in this section we 
have adopted the same definitions of normal galaxy and ULIRG given 
in their Table 2. 



Fig.|2j3 shows an updated version of the Lis.lL'^ cm\-n\ l umi- 
nosity ratio plot as function of Lir, where the data of'GS04b'have 
been replaced by the new HCN(l-O) data obtained with the 30m 
telescope for 16 LIRGs and ULIRGs. We also include in the plot 
the new HCN( l-O) dat a of the ULIRG IRAS 12112+0305 that 
was not part of lGS04bl 's sample. Quite noticeably, the inclusion 
of the new data has a significant effect on creating a tantalizing 
trend of LiR/LJ^^f^^j ratio with Lir. This trend, hardly evident 
in Fig. 12^, can be accommodated by a power law with indexes 
Of = 0.21 + 0.05 or a = 0.26 + 0.05, provided that the limits 
on ^HCNci-o) considered or discarded in the fit, respectively. 
With these new values, we derive a significant increase in the 
^iR / ^HCN( 1 -0) luminosity ratio that now goes from -900 Lq L' " ' 
for normal galaxies to ~17OOL0L' ' for ULIRGs. This trans- 
lates into a ~90% increase in Lir/L[j^j^^j over ~2.5 orders of 
magnitude in Lir, i.e., a factor of 4 larger than the corresponding 
increase derived from^S04b's data. 

The reasons explaining the changing picture that results from 
the comparison of Fig. |2^ and Fig. |2j5, are found on a sim- 
ple observational fact. A high percentage (-60%) of the ob- 
jects newly observed with the 30m telescope show HCN fluxes 
which are, on average, a factor of 2 lower than those reported by 
GS04b, Objects like Arp 193 are paradigmatic in this respect: 
the HCN(l-O) flux of Arp 193 in Table[f is a factor of 5 lower 
than that reported by lGS04bl Moreover, in other LIRGs and 
ULIRGs (like NGC 6240, Mrk 231, Mrk 273 and IRAS 17208- 
0014), the new HCN(l-O) single-dish fluxes reported in this pa- 
per are in much better agreement with the values derived from 
interferometer maps of the same sources (e.g., iTacconi et al.l 
IT99I lNakanishietani200l llmanishiet al]l2006l l2007h . Such 
an agreement is expected in the case of ULIRGs as the bulk 
of their molecular gas traced by CO lines and, very likely, of 
their dense molecular gas content traced by the HCN(l-O) line, 
is known to be concentra ted typically in their central kpc region 
(iDownes & Solomonll998i) . In this scenario the chances that in- 
terferometer maps may filter out a high percentage of the total 
HCN(l-O) emission in ULIRGs are scarce. Our results comfort- 
ingly fit this picture. 
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Fig. 2. (a) Lir/LJj(^(^j[ luminosity ratio as a function of Lir in the sample of normal galaxies, LIRGs and ULIRGs of lGS04bl 
Arrows indicate upper and lower limits to i'ecNci-O)- solid (dotted) line indicates the orthogonal regression fit calculated for the 
full sample of galaxies, excluding (including) upper and lower h mits. (b) Same as (a), but now the filled circles identify the new 
HCN(l-O) observations of LIRGs and ULIRGs i-eplacing |GS04bf s old data points. 



Many of t he conflicting s o urces were obtained with the 30m 
telescope by ISolomon et alJ (Il992h in several observing runs 
during a period going from 1988 to 1991, and were later com- 
piled in GS04b. In this context, it is worth noting that the cal- 
ibration scale of the 30m telescope is more accurate nowadays 
than two decades ago, owing to the new calibration unit and to 
the strong rejection of the image side-band of the new genera- 
tion of receivers installed in the late 1990s (IRAM Newsletter 
No. 38). In particular, the image side-band rejection used to be 
less accurately determined around 1990 that it is now, and fail- 
ure to set the proper sideband rejection could easily overesti- 
mate the calibration by up to a factor of 2 (cf. e.g. Kramer 1997, 
IRAM Report on Calibration of spectral line data). Furthermore, 
a significant improvement in the quality of the surface of the 
30m telescope, due to its re-adjustment in 1997 after holography 
measurements were carried out, has provided a better knowl- 
edge of the efficiencies Fgff and Bgff both governing the cali- 
bration scale of the data, and that are automatically set in the an- 
tenna control program since 2000 (IRAM Annual Report 2000). 
Finally, the pointing accuracy of the telescope is also much bet- 
ter today, due to a better knowledge, monitoring and control of 
the thermal and mechanical behavior of the antenna, now both 
taken into account in the pointing model (IRAM Annual Report 
2002). Altogether, we co nclude that the disagreement between 
the HCN fluxes given by ISolomon et al.l fl99Z) and those used 
in this work can be very likely attributed to occasional calibra- 
tion problems in the old data. 

In the following (Sect. |4|i we explore to what extent the ob- 
served trend of Lir/L^^j^^j with Lir can be interpreted as the 
signature of a change in the star formation laws in galaxies. 



4. Star formation efficiency and Kennicutt-Schmidt 
laws 

4.1. Working hypotheses 

In order to study how the SFR scales with the properties of the 
dense gas it is necessary to build a large sample of galaxies 
where Sdense and Ssfr can be both estimated. The objective is 




Fig. 3. Lir/Lfir versus Lft r for the sample of normal galax- 
ies, LIRGs and ULIRGs of lGS04bl (open circles), our sample 
of LIRGs and ULIRGs (filled circles) a nd a sample of inf rared- 
excess Palomar-Green QSOs (crosses: lEvans et aPlZOOd) . The 
solid line indicates the average Lir/Lfir luminosity ratio for nor- 
mal galaxies (Lfir < IO" Lq); the dotted lines show the range of 
the + standard deviation from the mean (~1.3+0.15). Infrared- 
excess Palomar-Green QSOs show the highest Lir/Lfir ratios 
due to the contribution of the AGN to the total MIR emission of 
the galaxy. 



to derive an updated version of the KS law for the dense gas us- 
ing this sample (Sect. |4.4| i. As a first approach to the problem, 
we make two assumptions: 

- First, we consider that L^j^j^^j is a good quantitative tracer 
of the total dense molecular gas mass (Mdense) of a galaxy in- 
dependently of Lfir, i-c, we implicitly assume that L^j^j^^j 
can be translated into Mdense using a universal conversion 
factor. However, the validity of this assumption is revised in 
Sect. |5]in the light of the new results presented in this work. 
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LpiR [ L©] LpiR [ Lq] 

Fig. 4. (a) Lfir/^^hcnci-O) luminosity ratio as a function of Lfir for the full sample of galaxies described in Sec. 14.21 Symbols are 
as in Fig. [3] The complete sample of high-z galaxies with available HCN observations is indicated with grey triangles dGao et al.l 
l2007i and references therein). Arrows indicate upper and lower limits to iecNci-O)- Horizontal lines mark the mean LpiR/incNci-O) 
luminosity ratios determined for normal galaxies, IR luminous galaxies and high-z objects, (b) Same as (a), but limits have not 
been represented. The solid (dotted) line indicates the orthogonal regression fit calculated for the full sample of galaxies, including 
(excluding) high-z objects. 



pointing to a significant overabundance of HCN in extreme 
luminous IR galaxies. 
- As a second assumption, we also adopt a universal conver- 
sion factor to derive the SFR o f a galaxy from its IR lu- 
minosity (e.g.. lKennicutll[r998h . However, this hypothesis 
can also be questioned if an AGN significantly contributes 
to the total IR power of a galaxy (LiR(8-1000;um)). In par- 
ticul ar, the AGN yield can be fairly large in the MIR range 
(e.g.. lRowan-Robinson & Crawfordll989l) . This casts doubts 
on the results obtained in Sect. [3] where the reported trend 
of Lir/LJ^pj^jj with LiR may simply reflect an increasing 
AGN contribution to Lir, rather than a change in the star for- 
mation efficiency. 



In order to minimize the bias affecting the IR luminosity-to- 
SFR conversion, we use LpiR (40-500 fim) to estimate a bolomet- 
ric IR luminosity which is corrected at best from the AGN con- 
tribution. This IR luminosity (denoted if™) is used to estimate 
the SFR (Sect. |4.4| l. Lfir is first calculated for each galaxy of 
the sample fitting its rest frame MIR, FIR and sub-mm spectral 
energy distribution (SED) to a two grey body model. We then 
integrate the fit only over the 40-500 yum range. This method, 
that fits the whole IR SED of the galaxy, is more accurate than 
the standard Lfir determination, that only fits the two IRA S data 
points at 60jum and 100 yum (e.g. San ders & Mirabell 19961) . Lfir 
values are used in Sect. 14.31 to explore the Lhir-L^^j^i-j cor- 
relation in our sample. We then use Lfir (40-500 //m) to derive 
L^j^, re-scaling Lfir (40-5 00 /im) by a factor ~1.3. This scaling 
factor has been determined from the LiR(8-1000jum)/LFiR(40- 
500 i^ m) avera ge value derived for the sample of normal galax- 
ies of lGS04bl The underlying assumption in this calculation is 
that the AGN contribution is far less substantial in normal galax- 
ies than in ULIRGs or quasars (see Fig. [3]). Therefore we can 
better estimate the bolometric correction that should be applied 
to derive Lj™ from Lfir (40-5 00 jum) using normal galaxies as 
templates. 



4.2. Sample 

We have compiled a sample of 88 galaxies with published FIR 
and HCN observations in the literature. This compilation in- 
cludes the 17 LIRGs and ULIRGs with new HCN(l-O) obser- 
vations presented in this pap er and t he sample of normal galax- 
ies, LIRGs and ULIRGs of GS04b'. We also include a sample 
of infrared-excess Palomar-Green QSOs (Evans et al. 2006) and 
the comple te sample of hig h-z galaxies with available HCN ob- 
servations dGao et al.ll2007l and references therein). 

Because of their redshift, some of these galaxies have not 
been observed in HCN(l-O), but in HCN(2-1). For these sources 
^HCN(i-{)) computed assuming that their rotational line lu- 
minosity ratios are L'y^c. N(2-u/^'ucN( \ -m ~ similar to the 
mean value measured by lKrips et al.l (12007 ) in a sample of local 
galaxies, including a few LIRGs and ULIRGs. For complete- 
ness, we have also included in the sample the z - 4 quasar 
APM 08279+5255, observed in the J=5^ rotational transition 
of HCN (Wag g et al. 2005). To derive its HCN(l-O) luminos- 
ity, we adopt L'^c N(^~'^^/ ^'rcn(]^o\ ~ Q- ^' assuming the physical 
conditions used in lGarcfa-Burillo et al.l (l2006h for this source. 



4.3. The star formation efficiency and tiie LnR-L[^^~^^^^ 
correiation 

We represent in Fig.|4]the Lfir/LJ^^j^j^j luminosity ratio versus 
Lfir for the full sample of galaxies defined in Sect. 14.21 Under 
the assumptions made in Sect. 14. II the Lpir/L'-^^^^^^ q^ luminosity 
ratio can be interpreted as a measure of the star formation effi- 
ciency of the dense gas (SFEdense - SFR/Mdense)- In Fig.|4^ we 
can see that SFEdense increases with Lfir from normal galaxies to 
LIRGs & ULIRGs and high-z objects. This result confirms that 
SFEdense is higher in high-z galaxies compared to normal galax- 
ies (Gracia-Carpio et al. 2007; Gao et al. 2007; Riechers et al] 
120071) and clearly extends this trend for the first time to the lumi- 
nosity range of local universe LIRGs and ULIRGs. An orthogo- 
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nal regression fit to the full sample of galaxies, excluding (lower 
or upper) limits, gives (solid line in Fig.|4j3): 



log 



'HCN(I-0)y 



= (0.24 + 0.04) logLpiR + (0.17 + 0.42) (1) 



or 



-^FIR 



L' 



1.5 L 



0.24 
FIR 



(2) 



HCN(l-O) 



We note that a similar regression fit is obtained if high-z 
galaxies are not included (dotted line). The latter implies that 
the regression fit parameters are not at all constrained by high-z 
galaxies. 

Another illustration of this result is shown in Fig. |5] where 
we have plotted the LFiR-i'HCN(i-O) correlation for the full sample 
of galaxies. An orthogonal regression fit, excluding limits, gives 
(solid line in Fig. |5]i: 

logLpiR = (1.23 + 0.06) log L^jcNd-o) + (0-97 + 0.46) (3) 

or Lfir - 9 L'hcn(i-O) (4) 

Similarly to the previous case, excluding high-z galaxies 
from the fit (dotted line) has virtually no effect on the result. 
In summary, the best-fit correlation found between LpiR and 
L' 



HCN(l-()) 



is seen to be superlinear, contrary to previous claims 
in the literature (lGS04atlGS04bl) . The derived power index is in 
either case significantly larger than unity: ~1.2. 



4.4. The Kennicutt-Schmidt law 

In order to test if a similar superlinear behavior holds for the 
KS-law of the dense gas, we have determined Sdense and Ssfr, 
using the HCN(l-O) and IR luminosities derived above, and the 
molecular gas size estimates from published CO interferometer 
maps, that are available for most of the sources. 

We first t r anslat e Lfl^ into SFR using the same factor used 
bv lKennicuttl(ll998lPI : 



SFR[Moyr '] 



1.7xlO-'°Lf™[Lo] 



(5) 



Similarly, we derive Mdense from ^^hcnci-O) assuming the con- 
version factor of lGS04al : 



Mdense [Mq] ^ 10 L[ 



'HCN(l-O) 



[Kkms^'pc^] 



(6) 



Finally, to obtain Ssfr and Sdense we normalize SFR and 
Mdense using the size of the molecular gas distribution of the 
sources. In Fig. |6^ we have represented the derived Ssfr and 
^dense for the full Sample of galaxies. Ssfr and Sdense follow 
a very tight correlation over more than 4 orders of magnitude 
in Xdense- An orthogonal regression fit to the data results in a 
KS-law of the dense molecular gas with a power index - 
1.12 + 0.04: 



10^' 



10 



10 



10' 



10^' 




1.23 + 0.06 
FIR °= l-HCN 



10^ 



L'HCN(i-o)[Kkm s"^ pc^] 



Fig. 5. Lfir versus Lhcn(i-O) using the full sample of galaxies. 
Arrows indicate upper and lower limits to i'HCN(i-o)- solid 
(dotted) line represents the orthogonal regression fit to the full 
sample, including (excluding) high-z objects. As in Fig.|4] limits 
are not taken into account in the regression fits. Symbols are as 
in Figs. |3] and m 



The fit of the KS-law to the whole sample is still significantly 
nonlinear (A^ > 1). A close examination to the residuals from this 
global fit indicates that the population of normal galaxies would 
be better described by a shallower index power law. As illus- 
trated in the lower panel of Fig. [6^, residuals for galaxies with 
Lfir < 10'' Lq show a systematic trend with Edense from 0.5 dex 
at 10 M0 pc"^ to -0.5 dex at 10^ Mq pc"^. To investigate if there 
is a turn upward in the KS law at Lfir ^ 10" Lq, we split the 
sample into normal (Lfir < lO" Lq) and IR luminous galaxies 
(Lfir > 10" Lq) and made a two-function power law fit to the 
whole sample. The result of this fit, displayed in Figs. |6]3 and 
|6};, indicates that, although the two galaxy populations follow 
very well defined KS power laws over 3 decades in Zsfr and 
2^dense, their characteristic power indexes are substantially differ- 
ent: = 0.88 + 0.07 for normal galaxies, while = 1.14 + 0.06 
for IR luminous galaxies. The extrapolation to higher Edense of 
the KS-law fitting normal galaxies falls short of explaining the 
typical values of Esfr in IR luminous galaxies by up to a factor 
~6-7 (see Fig.|6};). 

In order to evaluate the significance of the two-function 
power law fit versus that of a single power law, values 
have been derived, resulting in an improvement when the two- 
function power law fit is adopted: the value decreases from 
1.9 to 1.4. We thus propose that the relation between star for- 
mation and the properties of the dense molecular gas in galaxies 
can be better described by two KS-laws: 



logZsFR = (1.12 + 0.04)logZde„se + (-2.10 + 0.12) (7) orEsFR-0.02Z°fl 



logEsFR = (0.88 + 0.07) log Edense + (-1-68 +0.15) 



or EsFR- 0.008 E^,'l 



(8) 



^ Note that lKennicutll ( |1998l) assumed that the IR luminosity appear- 
ing on the right side of his equation, originally Lir(8-1000 yum), should 
correspond ideally to radiation coming from star formation reprocessed 
by dust; as argued in Sect. 14. II Lj™ values derived in the present work 
are a more educated guess of the SFR than Lir(8-1000 yum), however. 



for galaxies with Lfir < 10" Lq, and: 

logEsFR = (1.14 + 0.06) logSdense + (-2.12 +0.21) 

or EsFR- 0.008 



(9) 
(10) 



(11) 



(12) 



for local and high-z IR luminous galaxies with Lfir ^ 10" L, 
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Fig. 6. (a) Surface density of star formation rate, Zsfr, against surface density of the dense molecular gas mass as traced by the 
HCN(l-O) line, Xdense, for normal (open squares) and IR luminous galaxies (filled squares) for which the molecular gas size is 
available. Arrows indicate upper and lower limits to Xdense- The solid line is the orthogonal regression fit to the full sample of 
galaxies. In (b) and (c) we separate between normal and IR luminous galaxies, respectively. The dashed line is the orthogonal 
regression fit to galaxies with LpiR < lO^'L© ^^'^ the dotted line is a similar fit to galaxies with LpiR > IO^Lq' including one 
Palomar-Green QSO and five high-z galaxies. Upper and lower limits are not considered in the regression fits. Lower panels show 
the residuals of the fits in logarithmical units. 



In the previous analysis the CO distribution size is used to 
evaluate the surface densities, because it is available for far more 
objects than the dense gas distribution sampled with HCN. Any 
bias due to an overestimation of the molecular gas size corre- 
sponding to the dense gas, here taken from CO, should have a 
negligible effect on the reported tilt of the regression fits for the 
two families of galaxies. The underlying reason is that both pa- 
rameters, SsFR and Xdense, are affected to the same extent by this 
bias. Should a residual effect exist, correcting for it would even 
reinforce the magnitude of the tilt, as in low luminosity objects 
the CO distribution tends to be more extended than the dense 
gas distribution whereas for high luminosity objects, although 
limited to some cases, there is observational evidence that both 
distributions are more similar in extent. 

In the following sections we discuss the potential bias that 
we may have introduced in estimating Ssfr and Zjense from LpiR 
and i'HCN(i-O) according to the assumptions of Sect. 14. 11 and how 
they may have affected the results presented in this section. 



5. Excitation and chemistry of dense molecular gas 
in IR luminous galaxies 

IGC0 6' analyzed the data issued from the first HCO^ survey 
of LIRGs and ULIRGs, finding significant evidence that the 
HCN(1-0)/HCO^(1-0) luminosity ratio increases with LpiR in 
infrare d luminous galaxies (Lfir > 10" L©; see Fig. 2a of 
IGC06I ' s paper). Fig. shows an updated version of the same 
plot based on the new HCN(l-O) data presented in this paper. 
The inclusion of this new data set corroborates the existence of 
a statistically significant trend of the HCN(1-0)/HCO+(1-0) lu- 
minosity ratio with Lfir for Lfir > 10" Lq. As argued in Sect.|3] 



the greatly improved quality of the new HCN data used to build 
up Fig. gives a higher significance to this result compared 
to our previous finding (GC06). Of particular note is the recent 
i ndepe ndent confirmation of a similar result by Jmanishi et al.l 
(l2007h . who studied a sample of infrared luminous galaxies us- 
ing the Nobeyama millimeter array. 

Taken at face value, the reality of the trend shown in Fig.|7t 
suggests that the overall excitation and/or the chemical proper- 
ties of the dense molecular gas change on average with Lfir. 
The inclusion into this analysis of the new J=3-2 line data of 
HCN and HCO^, obtained for 10 galaxies in our sample, can fur- 
ther constrain the physical parameters (ne, and Lk) and chemical 
abundance ratios ([HCN]/[HCO^]) describing the average prop- 
erties of the dense molecular gas in infrared luminous galaxies. 
This characterization can be achieved through a radiative trans- 
fer modelling of the line ratios. 

5.1. One-phase Large Velocity Gradient (LVG) models 

Without further observational constraints at hand, the simplest 
approach consists of using a one-phase Large Velocity Gradient 
(LVG) scheme to fit the three independent intensity ratios (R,) 
which are derived from the line intensities measured in our sur- 
vey. To evaluate the goodness of the fits we have applied a stan- 
dard ;^'^-test to the obtained solutions. We estimate the;^f^ of the 
3-parameter fit from the expression: 

, ^ {Ri - ™d'=l)2 

Z ^ 2 (13) 

where R™°^'^'^ are the best-fit model predictions for and cr, 
are the estimated uncertainties on R,\ values qI < li help to 
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Source 


HCN{l-0) 
HCO+(1-0) 


HCN(3-2) 
HC0+(3-2) 


HC0+(3-2) 
HCO+(l-0) 


IRAS 17208-0014 


1.49 


±0.19 


1.56 ± 0.22 


0.37 ± 0.06 


Mrk231 


1.05 


±0.11 


0.97 ± 0.16 


0.27 ± 0.04 


IRAS 12112+0305 


1.70 


±0.53 






Arp 220 


2.19 


±0.13 


4.00 ± 0.23 


0.22 ± 0.02 


Mrk 273 


0.98 


±0.15 


1.07 ± 0.23 


0.45 ± 0.06 


IRAS 23365+3604 


1.52 


±0.44 






UGC 05101 


1.89 


±0.38 




0.56 ± 0.12 


VllZw31 


0.87 


± 0.18 




0.34 ± 0.10 


NGC 6240 


0.54 


±0.05 


2.19 ± 0.31 


0.20 ± 0.03 


Arp 35 


1.09 


± 0.17 




n on 1 n no 
U./9 ± U.Uo 


Arp 193 


0.63 


±0.07 


0.37 ± 0.07 


0.34 ± 0.04 


NGC 695 


0.70 


±0.16 




<0.29 


Arp 299 A 


0.53 


±0.04 


0.33 ± 0.09 


0.28 ± 0.03 


Arp299B+C 


0.62 


±0.06 


<0.52 


0.21 ± 0.03 


NGC 7469 


0.68 


± 0.04 


1.34 ±0.32 


0.16 ±0.03 


Mrk 331 


0.79 


± 0.07 


0.54 ± 0.16 


0.26 ± 0.05 


NGC 7771 


1.04 


± 0.09 


1.66 ± 0.42 


0.11 ± 0.02 



identify good fits. We have chosen to fit the following ratios: 
HCN(1-0)/HCO+(1-0), HCN(3-2)/HCO+(3-2) and HCO+(3- 
2)/HCO^(l-0). These ratios, listed in Table [3] are displayed in 
Figs. |7^-to-c as a function of Lfir for the infrared luminous 
galaxies of our sample. Being a one-phase scheme, we have to 
assume that the beam filling factor of the emission is common 
for all lines. The LVG filling factors implied by best-fit solutions 
have been confronted on a case-by-case basis with the values es- 
timated from interferometer CO maps that exist for all the galax- 
ies of our sample (TablefTli. This helps to discard LVG solutions 
implying beam filling factors (/yfin) for HCN and HCO^ lines that 
exceed those derived from high spatial resolution CO maps. 

In summary, we find that, for any common set of physical 
parameters (hh, and Ty) simultaneously fitting all the line ratios, 
we require [HCN]/[HCO^] abundance ratios to be >5 for a sig- 
nificant number of LIRGs and ULIRGs in our sample (Table|4]l. 
This is illustrated in the upper and middle panels of Fig. |8] that 
show the LVG model predictions of the three considered line 
ratios as a function of 11^2, the HCO^ column density per ve- 
locity interval (A^hco+/AV) and the [HCN]/[HCO+] abundance 
ratio. Within the one-phase scenario, HCN/HCO^ luminosity 
ratios >1 in the J=l-0 and 3-2 lines cannot be explained with 
[HCN]/[HCO+] < 2 for the range of HCO+(3-2)/HCO+(l-0) 
luminosity ratios observed in our sample of LIRGs and ULIRGs, 
requiring higher abundances of HCN relative to HCO^ (see also 
Fig. |9j. This result is mostly independent of the value of Ty 
adopted in the models (values in the range Tk - 20 to lOOK 
have been explored in our search for the best fit solution) and it 
is a direct consequence of the critical densities of the rotational 
lines of HCN being a factor of ~6 higher than those of HCO^ . 
Quite interestingly, the case for overabundant HCN is more com- 
pelling at high Lfir- 

The bottom panels of Fig.[8]illustrate the best fit solutions ob- 
tained for Tk = 60K in three galaxies: Arp 299 A, IRAS 17208- 
0014 and Arp 220. These are fair representatives of the three dif- 
ferent categories of solutions identified in our sample of LIRGs 
and ULIRGs (Table Hand Fig.EJ. In Arp 220, analyzed in de- 
tail by Gracia-Carpio et al. (in preparation) with the additional 
input provided by the H'^CN and H'^CO^ rotational lines, it 
is estimated that [HCN]/[HCO+] > 30, in sharp contrast to 
the case of Arp 299 A where we estimate that [HCN]/[HCO+] 
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Fig. 7. (a) HCN(1-0)/HCO+(1-0) luminosity ratio as a function 
of Lfir in our sample of LIRGs and ULIRGs. The dashed line, 
that represents a linear regression fit to the data points, shows the 
trend, (b) and (c) Same as (a) but for the HCN(3-2)/HCO+(3- 
2) and HCO+(3-2)/HCO+(l-0) luminosity ratios, respectively. 
The positions of Arp 299 A (open triangle), IRAS 17208-0014 
(cross) and Arp 220 (open square) are highlighted. 



~ 2. IRAS 17208-0014 represents an intermediate case where 
[HCN]/[HCO+] ~ 10. 

As discussed in Sec. [1] different scenarios may account for 
HCN being overabundant in extreme ULIRGs. But more rele- 
vant to the discussion in this paper is the fact that high HCN 



abundances call for a lower L' 



■to-Mdense convcrsion fac- 
tor (Xhcn), particularly at high Lfir. If the emission of HCN(1- 



HCN(I-{)) 
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Fig. 8. Results of the one-phase LVG calculations for the HCN(1-0)/HCO+(1-0) (black contours, upper panels) and HCN(3- 
2)/HCO^(3-2) line ratios (black contours, middle panels) as a function of the HCO^ column density per velocity interval 
(A^HCO+/Ay), the H2 volume density («h,) and the [HCN]/[HCO+] relative abundance {- 2, 10 and 30 from left to right). The 
grey area indicates the range of HCO^(3-2)/HCO^(l-0) line ratios (white contours, upper and middle panels) observed in our 
sample of LlRGs and ULIRGs. In the bottom panels we show the one-phase LVG best fit solution of three representative sources: 
Arp 299 A, IRAS 17208-0014 and Arp 220. Filled, dashed, and dot-dashed contours represent, respectively, the regions where the 
predicted HCO+(3-2)/HCO+(l-0), HCN(1-0)/HCO+(1-0) and HCN(3-2)/HCO+(3-2) Hne ratios are equal to the observational 
values given in Table |3] Thinner lines indicate one-cr deviations from these ratios. Grey areas correspond to regions where < 3, 
1.5 and 0.5, from light grey to black. The kinetic temperature was fixed to 60K in all the calculations. 

0) is optically thin (a very unlikely scenario), the value of Xhcn tively evaluated from the LVG results reported above. If we as- 
should scale down as 1/[HCN]. It must be emphasized that sume that the HCN(l-O) line comes from an ensemble of viri- 
even in the case of optically thick lines (a more realistic see- alized, non-shadowing molecular clumps with optically thick 
nario, favored by our LVG analysis), a significant change in 
^HCN, due to the different properties of the emitti ng gas, can- 
not be excluded either (see also discussion in Papadopoulos] 
12007; .Papadopoulos et al.,,2007,) . This change can be quantita- 
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[HCN]/[HC01 

Fig. 9. Results of the one-phase LVG calculations for 
the HCN(1-0)/HCO+(1-0) (black contours) and HCN(3- 
2)/HCO^(3-2) line ratios (grey contours) as a function of 
the [HCN]/[HCO^] abundance ratio and the H2 volume den- 
sity (hhJ. The HCO^ column density per velocity interval 
(A^HCO+/Ay) and the kinetic temperature (Tk) have been fixed 
to 10'^ cm"- km"' s and 60 K, respectively, to illustrate the 
range of solutions for [HCN]/[HCO^]. Symbols indicate the 
best fit solutions of the three representative sources discussed in 
Sect.EU Arp 299 A (open triangle), IRAS 17208-0014 (cross) 
and Arp 220 (open square). 



emissio n, Xhcn can be estima ted from the expression first de- 
rived bv lDickman et alj (Il986h . as: 



1/2 



HCN 



2.1 X 



7"hcn 



(14) 



where rtu^ and Thcn represent respectively the density and 
brightness temperature derived from the LVG solution for the 
HCN(l-O) line. Based on the output from the best fit models 
(summarized in Tabled, we estimate that the Xucn value is ~4.5 
times lower at LpiR ~ 10'^ L0 compared to the value typically 
found at LpiR ~ 10" Lq. While a better quantitative assessment 
of how Xhcn changes on average with LpiR would require multi- 
phase radiative transfer models, a direct implication of this find- 
ing is that the reported trend in the Lfir/^hcn ^''^^^'^ a function 
of LpiR, discussed in Sect. 21 may be hiding a potentially more 
dramatic change of the SFEdense- 



5.2. Two-phase Large Velocity Gradient (LVG) models 

It is noteworthy that the simplest one-phase LVG approach de- 
veloped in Sect. 15. H is akeady able to fit the three line ratios in 
most of the galaxies of our sample with a satisfactorily low 
figure (<3), except for three sources: in NGC 6240, NGC 7469 
and NGC 7771, the derived;^- is too large (3>3). More sophisti- 
cated models, accounting for the multi-phase nature of molecu- 
lar gas, can be used to try to significantly improve the fit to the 
three molecular line ratios observed in the galaxies of our sample 
in general, and in the three conflicting sources referred to above 
in particular In addition, with the help of multi-phase models, 
we can explore if the main conclusions of Sect. 15.11 favoring a 
significant change of Xhcn with LpiR are a mere artefact of one- 
phase LVG schemes. One-phase models may not properly take 



into account that the critical densities of the rotational lines of 
HCN and HCO^ for the same J level differ by a factor ~6. The 
latter would favor that the emission of these two lines may come, 
on average, from two distinct phases characterized by different 
densities. 

The main drawback of multi-phase models is that the num- 
ber of free parameters to explore a priori is exceedingly large, 
however The solutions are usually degenerate and the relevant 
parameters cannot be well constrained if the available number of 
data points to fit is reduced compared to the degrees of freedom. 
In our case, the constraints are given by the three line ratios in- 
dependently measured in each galaxy, and also, by the the upper 
limit on the size of the source derived from CO interferometer 
maps (9co in Table [TJ. 

In order to minimize the number of free parameters, we have 
reduced the number of phases to two: a hot (T^ = 80 K) and 
dense («h, = 10'' cm"-') phase, which would be more directly 
associated with massive star formation, and a cold (Tk - 25 K) 
and comparatively more diffuse phase (nu^ = 10"^ cm"^), which 
would account for a quiescent reservoir of molecular gas. This 
dichotomy in reflects the two-temperature fit of the aver- 
age IR SED of the galaxies in our sample. Furthermore, the 
density range explored by the iien.se and diffuse phases pur- 
posely encompasses the values of the effective critical densities 
of HCN( 1-0) and HCO+(1-0) respectively. The [HCN] / [HCO+] 
abundance ratios are also fixed for the two phases as follows. We 
adopt a canonical value of [HCN]/[HCO^] - 1 for the diffuse 
phase, typical of quiescent molecular dark clouds in our Ga laxy 
( lOhishi et alJll992l: iPratap et al.lll997l: fPickens et alj|2000l) . In 
contrast, we fix [HCN]/[HCO^] - 10^ in the dense phase, a 
value typical of galactic sta r forming regions where hot-core 
like c hemistry can develop jjohansson et al.|[l984t iBlake et aTl 
|1987| ). This choice of physical and chemical parameters for the 
two phases is expected to alleviate the need of enhancing the 
HCN global abundanc^ with respect to that of HCO^ in some 
of the sources, a result found in Sect. 15. II The robustness of this 
result and its implication for the constancy of Xhcn can thus be 
tested. A side effect of purposely fixing some of the parameters 
of the fit, for the reasons explained above, is that the resulting 
may be higher, i.e, formally worse, in the two-phase fit. 

In the fitting procedure, the relative filling factor of the two 
phases (/ for the dense gas and 1-/ for the diffuse phase) is 
allowed to vary. We also allow to freely fit the column density 
of HCN in the two phases with the constrain that [HCN] > 
[HCN]'''"^"^'^, in order to be consistent with the values derived in 
our Galaxy for cold dark clouds and star forming regions. As 
in Sect. 15.11 we only keep the LVG solutions with rjan values 
for HCN and HCO^ lines that are below those derived for CO. 
Table|4]summarizes the results of the two-phase model fitting. To 
ease the comparison with one-phase models, we have derived the 
global [HCN]/[HCO^] abundance ratios and Xucn conversion 
factors characterizing the best fit two-phase solutions. Global es- 
timates are obtained by averaging the contributions from the two 
phases, as follows: 



[HCN] fN\ 



f\ A7diffuse 



[HCO+] fN^co* + (1 



f\ AT diffuse 



ry^dense vdense _i_ / I _ ^\T^diffuse ^diffuse 



"^HCN 



XNT^ diffuse 
■ HCN 



(15) 



(16) 



^ We mean by global HCN abundance, the average value derived 
from the two-phase model. 
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Table 4. One-phase and two-phase LVG best fit solutions for the sample of LIRGs and ULIRGs detected in the J=l-0 and 3-2 
lines of HCN and HCO^. values of the fits are given in Col. 2. The predicted line ratios are in Cols. 3 to 5, followed by the 



[HCN]/[HCO^] abundance ratio (Col. 6), the L^j^j^^j Q^-to-Mdense conversion factor (Col. 7) and the estimated molecular source size 
at FWHM (Col. 8). For the two-phase LVG model, the relative filling factor of the dense phase is also given (Col. 9). Galaxies are 
ordered with decreasing Lfir. L' - Kkm s ' pc^. ^ Global values computed according to Equations[T5]and[T6l 

One-phase LVG model 



Source 



HCN(l-O) 
HCO+(I-0) 



HCN(3-2) 
HC0+(3-2) 



HC0"^(3-2) 
HCO+(1-0) 



[HCN] 
[HCO+] 



^HCN 

(MqL'- 



0s 



IRAS 17208-0014 


0.0 


1.49 


1.58 


0.37 


10.5 


45 


0.67 


Mrlc231 


2.7 


1.04 


1.08 


0.33 


6.8 


26 


0.85 


Arp220 


1.3 


2.32 


3.92 


0.23 


28.0 


22 


1.24 


Mrlc273 


0.0 


0.99 


1.04 


0.45 


6.6 


17 


0.39 


NGC 6240 


33.3 


0.55 


0.38 


0.20 


2.5 


127 


1.69 


Arp 193 


0.0 


0.63 


0.37 


0.34 


2.2 


130 


1.14 


Arp299A 


0.0 


0.53 


0.33 


0.28 


2.0 


127 


1.56 


NGC 7469 


6.5 


0.69 


0.54 


0.17 


3.5 


147 


2.09 


Mrk331 


0.0 


0.79 


0.55 


0.26 


3.4 


123 


1.22 


NGC 7771 


12.1 


1.01 


0.85 


0.18 


5.4 


129 


2.12 



Two phase LVG model 



Source 




HCN(l-O) 
HCO+(1-0) 


HCN(3-2) 
HCO+(3-2) 


HCO+(3-2) 
HCO+(1-0) 


[HCN] a 
[HC0+] 


^HCn'' fs 

(MqL'-') (") 


/ 

(%) 


IRAS 17208-0014 


4.8 


1.16 


1.79 


0.42 


3.1 


23 


0.60 


17.4 


Mrk231 


11.8 


0.84 


1.24 


0.37 


2.2 


31 


0.85 


10.5 


Arp 220 


17.9 


1.77 


4.46 


0.25 


4.2 


35 


1.28 


24.2 


Mrk 273 


0.4 


0.90 


1.15 


0.46 


2.4 


23 


0.42 


12.0 


NGC 6240 


0.9 


0.57 


1.97 


0.19 


49.6 


132 


1.73 


1.6 


Arp 193 


4.9 


0.48 


0.43 


0.36 


1.3 


44 


0.94 


2.6 


Arp 299 A 


7.7 


0.47 


0.47 


0.31 


1.3 


56 


1.27 


2.7 


NGC 7469 


1.2 


0.67 


1.61 


0.18 


1.7 


88 


1.53 


7.0 


Mrk 331 


10.1 


0.64 


0.84 


0.34 


1.7 


40 


0.87 


6.3 


NGC 7771 


11.5 


0.85 


2.63 


0.15 


2.0 


103 


1.96 


9.4 



In the first equation, N^^^ and N^^^'^ represent the column 
densities of HCN derived for the dense and the diffuse phase 
model clouds, respectively. Similar symbols are used for HCO^ . 
In the second equation, ^hcn ^hcn factors 

calculated for the dense and the diffuse phase, respectively, ac- 
cording to Equation[T4] The brightness temperatures of HCN(1- 
0) in the dense (r^^N ) ^"^^ i" diffuse phases (r^^;^'') are 
used as weighting factors in each phase. 

First, as a bottom line conclusion, we see that the goodness 
of the fit of the two-phase models, measured by is worse 
compared to that achieved in the one-phase models. This applies 
to all sources, with the exception of NGC 6240 and NGC 7771 
(see Table m. Although with a large scatter, the best-fit models 
indicate that the filling factor (/) of the dense phase increases 
with Lfir. This result can be taken as an evidence that the aver- 
age density of molecular gas increases with Lfir. In addition, we 
find no significant difference in the [HCN]/[HCO^] abundance 
ratios, estimated by Equation[T5] between ULIRGs and LIRGs, 
i.e., we eliminate the need of enhancing the global abundance 
of HCN with respect to that of HCO^ at high Lfir, in contrast 
with the result found in Sect. 15.11 NGC 6240 is, on this respect, 
also an exception to the rule. The increasing trend of / with Lfir 
indicated by the best-fit solutions seems to neutralize the need 
of a similar trend of [HCN]/[HCO^] with Lfir in the two-phase 
model. However, we find that the conversion factor Xhcn, de- 
rived from Equation[T6l is lower at high Lfir, a conclusion simi- 
lar to that found in one-phase models. In particular, Xucn is ~2.5 



times lower at Lfir ~ 10'^ L0 compared to the value typically 
found at Lfir ~ 10" Lq. 

6. Conclusions 

In this paper we present observational evidence that the 
Lfir /Lhcn( 1-0) ratio, taken as proxy for SFEdense, is a factor 
~2-3 higher in galaxies categorized as IR luminous (Lfir > 
10" Lq) compared to normal galaxies. Local universe LIRGs 
and ULIRGs populate a region in the SFEdense diagram that 
lies between those occupied by normal and high-z IR luminous 
galaxies. The reported trend in the SFEjense derived from HCN 
data implies that there is a statistically significant turn upward in 
the KS law, Esfr ^ ^dense' ^^8^ ^fir: N changes from ~0.80- 
0.95 (for Lfir < 10" Lq) to -1.1-1.2 (for Lfir > 10" L©). 
Furthermore, our multiline analysis of HCN and HCO^ data in- 
dicates that XucN is ~3 times lower at high Lfir. This latter find- 
ing reinforces a scenario where the SFEdense niay be up to an 
order of magnitude higher in extreme IR luminous galaxies than 
in normal galaxies. 

Ba sed on the model developed by iKrumholz & McKe3 
(120051) . in which the SFR density (psfr) scales with the aver- 
age density of the gas as ~ p^'^f . Krumholz & ThompsonI (l2007h 
conclude that the power index of KS laws determined directly 
from observations would change depending on how the effective 
critical density of the tracer used to probe the star forming gas 
compares to the average density of the gas itself. In the particu- 
lar case of HCN(l-O) data, only in galaxies where pg^^j, exceeds a 
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few 10"^ cm"^ (i.e., the effective critical density of the HCN J=l- 
line), we would start recovering the expected superlinear be- 
havior of the universal KS law derived by Krumholz & McKee 
(l2005h . Our results, hinting at a break in the KS power law de- 
rived from HCN data around LpiR = 10" Lq, are in qualitative 
agreement with this picture. However, it remains to be proved 
if the one order of magnitude increase in the SFEdense over 3 
decades in LpiR derived from our data can be fitted by a univer- 
sa l KS law with a corresponding increase of Pg^j, (e.g., see Fig. 2 
of lKrumholz & Thompsonl2007h . Based on our two-phase LVG 
fits of LIRGs and ULIRGs molecular data, we do find tantaliz- 
ing evidence of higher p,^, at higher LpiR. Enlarging the number 
of detected molecules and transitions will be a key to better con- 
straining pg,,. 

Alternatively, the reported increase of the SFEdense at high 
LpiR may be due to the fact that star formation processes devi- 
ate from the KS recipes followed by normal galaxies. The high 
frequency of interactions can create a pressure enhanced envi- 
ronment in the ISM of LIRGs and ULIRGs that may specifically 
elevate SFEdense in these galaxies. On the other hand, an extreme 
ISM environment, characterized by high densities and tempera- 
tures, may also create a top-heavy Initial Stellar Mass Function 
(IMF) in luminous IR galaxies. This scenario has been suggested 
to ap ply to the nuclei of starbursts and to the center of the Galaxy 
(e.g..'Klessen et alll2007i) . With this assumption, the conversion 
factor between SFR and Lfir, mostly IMF-dependent, should be 
lowered at high Lfir, instead of being kept constant as supposed 
in Sect. |4] The SFEdense trend derived above may thus be reflect- 
ing a change of the IMF rather than a true variation of the SEE. 

Finally a significant contribution from an AGN source to 
Lpi R canno t be excluded, especially at the high luminosity end 
(cf. lGC06l) . This would produce a trend in the SFEdense sim- 
ilar to the one reported above but due to the AGN contam- 
ination being more severe at high Lh r- In this context, the 
discovery by iDownes & EckartI ilOO 't) of an embedded AGN 
source in the West nucleus of Arp 220 is noteworthy: in this 
ULIRG, considered to date as a prototype o f a starburst dom- 
inated IR galaxy (e.g., iGenzel et al] Il998l) . iDownes & EckartI 
( l2007h claim that -75% of flie total IR luminosity of the 
galaxy may come from the accretion disk of the AGN. The 
fact that we have also derived in Arp 220 the most extreme 
[HCN]/[HCO^] abundance ratio among all the galaxies of our 
sample, points to a causal link between HCN overabundance 
and AGN-driven chemistry. However, it is still a mostly con- 
troversial issue whether X-ray chemistry is able to enhance 
[HCN]/[HCO^] abundance ratios in AGNs to the le vel required 
by observations ( Usero et all 2004 : Meiierink& S paansI 120051: 
iGarcfa-Burillo et al.ll2006t iMeiierink et al.ii2007 ). Other related 
processes involving the evaporation of dust grain mantles and 
high temperature gas-phase reactions, likely efficient in the 
dense and hot molecular environments typical of extreme star- 
bursts, and probably also of AGNs, may be responsible of shap- 
ing^ the chemistry of the molecular gas in IR luminous galaxies 
(IGC06: Lahuis et al. 20C^. 

This paper shows how multiline analysis of molecular tracers 
of dense gas can contribute to study the star formation rate and 
star formation efficiency laws in starbursting galaxies, providing 
more elaborated tools to determine whether the SEE is enhanced 
or constant in the most extreme starbursts found in the local and 
the high-z universe. 
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